Bosonization approach to charge and spin dynamics of ID fermions 

with band-curvature 
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We consider one-dimensional (ID) spin-1/2 fermions in a clean quantum wire, with forward scat- 
tering interactions and a non- linear single-particle spectrum, ^fc = /2m where v is the Fermi 
velocity and 1/m is the band-curvature. We calculate the dynamical structure factor (DSF) of the 
model at small wave-vector q with the help of the bosonization technique. For spinless fermions, 
we show that, starting from the single-parametric spectrum: ui — u\q\, bosonization emulates the 
2-parametric excitation spectrum: uj = u\q\ ±g^/2m*, where m* decreases with increasing repulsive 
interactions. Moreover, away from the excitation-cone, i.e. w ^ u\q\, bosonization yields the 2-pair 
excitation continuum of the DSF. For spinful fermions, we show that the spin-charge coupling (SCC) 
due to band-curvature affects charge and spin DSF in an asymmetric way. For the charge DSF, SCC 
manifests as a two-peak structure: a charge peak at u = Up\q\ but also a spin peak at = y-cr\q\, as 
charge fluctuations may decay via chargeless spin-singlet excitations. For the magnetic DSF, SCC 
manifests as a continuous transfer of magnetic spectral weight to frequencies uj > ■Uo-|'/|, as spin 
fluctuations decay via pairs of chargeless spin and spinless charge-neutral excitations. 



I. INTRODUCTION 

Correlation functions of one-dimensional (ID) fermions 
in quantum wires are conveniently calculated within the 
Tomonaga-Luttinger (TL) model and with the help of 
the bosonization technique. The TL model assumes the 
linearity of the single-particle spectrum with respect to 
momentum: 

S,k — v{\k\ — kp) (Tomonaga — Luttinger), (1) 

where v is the Fermi velocity and kp the Fermi momen- 
tum, and interactions of the forward scattering type, so- 
called g2 and 54 processes in g-ology. The forward scat- 
tering nature of the interactions implies that the sys- 
tem has no spectral gap whereas the linearity of the 
spectrum yields Lorentz invariance. These assumptions 
make the model exactly soluble by mapping the interact- 
ing fermions onto free bosonic excitations or plasmons 
a consequence, e.g. transport, proper- 
ties of TL liquids are easily accessed with the help of 
bosonization, see the recent monographs [^l7|. 

Even though the TL model allows a standard descrip- 
tion of fermions in quantum wires, as Landau theory of 
Fermi liquids allows a standard description of 3D sys- 
tems, cf. the monograph 8], both assumptions of rel- 
ativistic single-particle spectrum and forward scatter- 
ing interactions often constitute an over simplified de- 
scription of ID fermionic systems. As a matter of fact, 
fermions hopping with an amplitude i on a ID lattice of 
parameter a have a single-particle spectrum given by: 



-2tcos{ka) (Hubbard). 



assume that the system has no spectral gap, e.g. is 
not in a Mott insulating phase. One may then focus 
on the low-energy properties of this model by expand- 
ing Eq. ^ around the Fermi points ±kp where the -I- 
sign refers to right-movers and the — sign to left-movers. 
In the lowest order, this yields the TL model, Eq. ([1]), 
with V — 2tasm{kpa). Because of the exact solubility of 
the TL model the low-energy properties of the Hubbard 
model are thus known exactly. On the other hand, ex- 
cept for the case of free fermions which is again exactly 
soluble, less is known about the high-energy properties of 
the Hubbard model. Considerable work has been devoted 
to exact studies related to the nature of its ground-state 
and its low-energy excitations. This has been done with 
the help of the Bethe Ansatz, cf. [t^lQ] for reviews, which 
allows a fully non-perturbative handling of the non-linear 
spectrum, Eq. and interaction effects. This technique 
is however extremely difficult to generalize to the calcu- 
lation of correlation functions. Formidable efforts are 
made to achieve this task but they are limited at present 
to the 1/2— filled case [ll| and to the XXZ Heisenberg 
spin-chain, see for a review and Ref. [2l[ for an ap- 
plication to the present problem (the XXZ spin-chain is 
equivalent to spinless fermions, cf. Ref. [l^ for a review 
on spin-chains and the Hubbard model). 

In this contribution we focus on a different approach 
to the high-energy properties of strongly correlated ID 
fermions. This amounts to start from the TL model, as- 
sume forward scattering interactions and introduce non- 
linear corrections to the linear single-particle spectrum. 
Restricting ourselves to the lowest order correction yields: 



(2) Ck=v{\k\-kp) + {\k\~kpf/2m + 0{{\k\-kp)^), (3) 



Adding interactions between such fermions leads to the 
ID Hubbard model, see the monograph ^9||, one of the 
most fundamental model of solid-state physics. In the 
absence of disorder and away from 1/2— filling we may 



where 1/m ~ 2ta^ cos{kpa) is the band-curvature or, 
equivalently, m is the band-mass. For Galilean invariant 
systems, there are no higher order corrections than the 
band-curvature. On the other hand, for lattice fermions. 
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Eq. ([3]) neglects an infinite number of irrelevant (in the 
renormalization group sense) terms. This means that, 
with respect to the Hubbard model Eq. ([2|), the model 
defined by Eq. ^ only takes into account the most rele- 
vant of these terms which, below half-filling, corresponds 
to the curvature of the band-spectrum. Such a simplified 
model then allows a non-perturbative treatment of band- 
curvature and to extract information about dynamical 
correlation functions of interacting ID fermions beyond 
the TL model and towards the Hubbard model. 

Historically, such a generalization of the Tomonaga 
model has first been considered by Michael Schick Q, 
see also (5lfl4[. More recently, this model has been con- 
sidered in relation with the damping of ID plasmons, 
Refs. fT5lM'20lf2ll, and the drag resistivity due to for- 
ward scattering between quantum wires, Ref. (IS] . Most 
of these studies focus on the density-density correlation 
function in the case of spinless fermions, see Refs. [l^, 
[i8,19;20;2L22,] or the equivalent XXZ Heisenberg spin- 
chain, see Refs. [20ll2l| . Various techniques have been 
considered to deal with this pro blem; mainly: fermionic 
fl6, 18, 19,20,24 field theory [U, Bethe Ansatz [U and 
numerical (Density-Matrix Renormalization Group) [2l| . 
Non-perturbative approaches to band-curvature may be 
found in [16^20,21]. In this paper we will focus on a 
field theory (bosonization) approach to the computation 
of density-density correlation functions and consider the 
case of fermions with spin. 

The paper is structured as follows. In Sec. [H] we will 
review some known results (for free as well as interacting 
20||2]| ) fermions and motivate the need for a bosoniza- 
tion approach. In Sec. IIIII the case of spinless fermions 
will be considered and the results compared with the ones 
found in the literature, cf. 20||2]| . In Sec. llVl the case of 
spin— 1/2 fermions will be considered following Ref. p^ . 
In Sec. |V]we will summarize our results and conclude. 



II. MOTIVATIONS FOR A BOSONIZATION 
APPROACH 




FIG. 1: Schematic view on the spectrum of excitations of 
ID spinless fermions. In the case of free fermions, a single- 
pair excitation continuum, S° , lies between the u)- and u!+ 
branches. In the presence of forward scattering interactions 
a multi-pair excitation continuum (S'^' at the 2— pair level) 
spreads above the uj+ branch and a divergent spectral weight 
appears along the cj- branch, cf. Fig[21 



2— parametric family appears: 



(5) 



In the (particle-hole) excitation spectrum of the system, 
this family bounds the region of decay of the coherent 
bosonic excitation, the plasmon which corresponds to a 
coherent particle-hole pair, into the incoherent contin- 
uum of single particle-hole pairs, see Fig. [T] This is 
best revealed by looking at the dynamic structure factor 
(DSF) of the fermions which corresponds to the dissipa- 
tive part of the retarded polarization operator: 



5(w,g) = -3H^(t^,<7). 
For free fermions, Eqs. ^ and (O yield (at T = 0): 



(6) 



\l\ 



e[ (w- 

q 



(7) 



The main object we will consider in the following is 
the density-density correlation function or polarization 
operator: 

H(x,r) = (p(x,t)p(0,0)), 

where p is the fermion density measured relative to a 
neutralizing jellium. 

For free-fermions this correlation function is known ex- 
actly with band-curvature. It is simply the polarization 
bubble, i.e. the convolution of two free fermionic Green's 
functions. In Fourier space it reads: 



HO(zc^,g) = — In 
T7q 



,0 \2 



(4) 



where the upper index refers to free-fermions and we have 
taken into account the spin degeneracy. In Eq. Q, the 



which, for a fixed momentum is non-zero and inde- 
pendent on frequency in the range: uj'^ < to < lu^. In 
particular, at the plasmon frequency, Eq. ([7]) yields a fi- 
nite result: 



S"{uj = vq,q) 



m 



(8) 



Moreover, the spectral width of the DSF yields the in- 
verse life-time of the bosonic excitations: 



m 



(9) 



On the other hand, expanding Eq. (jj]) in 1/to yields: 



U°iiu;,q) 
2 ^ 



(10) 



TTV ^-^ 2n - 

n=0 



(vq) 



2n+2 



[(iw)2 - (vq) 



2]2n+l ' 
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FIG. 2: Schematic view on the dynamical structure factor of 
interacting spinless fermions as a function of frequency uj for 
a given momentum q, e.g. the cut at g <C 2kF of Fig. [1] (at 
T = 0). The single-pair excitation continuum lies between 
the ui^ and lj+ branches. As known from Refs. [20ll2ll |. the 
spectral weight diverges at the lower-edge lu- and crosses over 
smoothly to a multi-pair continuum at the upper edge uj+. 



which reveals a basic difference with respect to the ex- 
act result of Eq. ([4]) : order by order the dissipative part 
is singular at the plasmon frequency: u = v\q\^ with 
5— function singularities. In particular, within the TL 
approximation (to — > 00): 



(11) 



which implies that plasmons have an infinite life-time, 
i.e. are free. The divergency of the DSF, in the limit 
TO — !■ 00, agrees with Eq. ([5]) but the nature of the di- 
vergency is different according to the order with which 
the limits to ^ 00 and uj ^ vq are taken. The fact that 
these limits do not commute witness the non-trivial role 
of the irrelevant band-curvature correction. Indeed, the 
latter cures the TL (5— function singularities, c/. Eq. ([5]), 
providing the exact line-shape of the free-fermion DSF, 
c/. Eq. ([7]), and the finite life-time of the elementary 
excitations, c/. Eq. Moreover, all these results are 
non-perturbative in band-curvature as they require sum- 
ming over an infinite number of singular terms. 

The case of interacting fermions has been the sub- 
ject of the recent literature. The results are summarized 
schematically on Fig. [21 The latter shows that the in- 
terplay between band-curvature and interactions for the 
DSF is two-fold. On the one hand, at high frequencies 
(w ^ vq)^ this interplay yields a multi-pair continuum to 
the DSF. The 2— pair continuum was first derived in Ref. 
[3 and subsequently in Refs. In Refs. [Hill 

fermionic techniques were used to derive the 2— pair con- 
tinuum whereas Ref. [2l| derived it with the help of 
the bosonization technique. This high frequency contin- 
uum may be accessed by perturbation theory in curvature 
and interactions. On the other hand, at low frequencies 
{uj « ug), a (power-law) divergency of the spectral weight 
appears at the lower-edge w_ and a smooth crossover to 
the multi-pair continuum takes place at the upper edge 
UJ+, [20,21] . Moreover, the width of the single-pair con- 
tinuum was shown to be affected by interactions [T6ll2]| : 



where to* gets smaller with increasing repulsive interac- 
tions, i.e. the plasmon peak becomes broader. Contrary 
to the high-frequency part the low-frequency regime re- 
quires a non-perturbative treatment of band-curvature. 
This was already the case for free fermions. Adding inter- 
actions of course considerably complicates the problem. 

IgI 2oI | by fermionic techniques 



This was achieved in Ref. 

and in Ref. ^2l| by Bethe Ansatz and numerical simula- 
tion. As can be seen from Fig.[2]the free-fermion box-like 
line-shape is therefore significantly modified by interac- 
tions, both at low and high frequencies. 

In this frame, our goal will be to calculate the DSF (at 
small q, i.e. the 2fci7'— part is not considered) with the 
help of the bosonization technique. We will start from a 
Hubbard chain below 1/2— filling and show in detail how 
the bosonization can be implemented when taking into 
account band-curvature corrections to the single-particle 
spectrum. Our first motivation for using bosonization 
comes from the fact that it allows to take forward scat- 
tering interactions among the fermions to all orders even 
though it is essentially perturbative in band-curvature. 
Following the recent Ref. [2l[ we will reconsider the case 
of spinless fermions. The exact results available in this 
case will allow us to discuss the efficiency and limita- 
tion of this technique towards a possible re-summation 
of band-curvature corrections. Our second motivation 
comes from the fact that bosonization allows to take 
quite straightforwardly into account the spin- 1/2 of the 
fermions. We will focus on the highly non-trivial coupling 
between spin and charge degrees of freedom due to band- 
curvature. Within bosonization, this was first noticed in 
Ref. |1J| that we shall follow (correlation functions were 



not the main issue of Ref. [lj|). A spin-charge coupling 



due to band-spectrum non-linearities was also reported 
in Ref. within a Bethe Ansatz approach to the ID 
Hubbard model (correlation functions are not yet com- 
putable with this technique). As we will detail below, 
field and lattice theories agree qualitatively. On this ba- 
sis, we will explore, with the help of the bosonization 
technique, the effect of this coupling on the charge and 
spin density correlation functions of the model. At this 
point, the reader which is interested more in our results 
than the derivations may refer directly to Sec. |Vl 



III. SPINLESS FERMIONS 
A. Model 

In this Section, we consider spinless fermions and give 
some details of the standard bosonization procedure. 
Starting with free fermions their Hamiltonian reads: 



H 



k 



5u) cx 



g 

m* 



(12) 



where c and are the annihilation and creation opera- 
tors for the fermions, is given by Eq. ([3]) and :: refers 
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to normal ordering. The bosonization procedure is im- 
plemented by first focusing on the low-energy sector near 
the Fermi points ±kp: 



+ikFX 



c{x) ~ il)+{x)e-'''^'' + i}}^{x)e 

where 7p± are slow chiral fields. The free fermion Hamil- 
tonian becomes: 



H — iv 



dx 



-— I dx 

2m 



.(13) 



The first term corresponds to the Dirac part related to 
the linear spectrum. The second term corresponds to the 
curvature correction. One may then go to the bosonic 
representation with the help of the following identity: 



ij±{x) 



1 



exp 



±i\/47r (p±{x) 



(14) 



relating the chiral fermionic fields, tp±, to the chiral 
bosonic fields, (p±. The latter may be expressed as a 
function of the total phase field and the momentum 



11;^ as follows: 



T / dx' U^{x') 

'J —OQ 



(15) 



The normal-ordering of vertex operators is conveniently 
taken into account with the help of a point-splitting tech- 
nique. This yields the following operator-product expan- 
sions (77 = ±): 



1 



V'^(a;)V'r,(a;) = -—^ d^Lp^ix), 



i>l{x)dxi>r,(x) = ' 

i>l{x)dxx'lpr,ix) = 



-i sgn(?7) : (da^iprjix))'^ 



(16) 



where the terms which vanish upon integration over space 
or by symmetry considerations have been neglected. The 
last line in Eq. displays the non-linearity of the 

bosonization arising from band-curvature. The fermionic 
Hamiltonian then takes over the following form: 



H 



V J dx [{d.,ip+f + id,^-f] 

2\/^ f r/a ^3 , /a ^3l 



3m 

with the densities: 

P+i^) = 7 



dx [{d^if+f + (d^cp^f] , (17) 



dxip+{x), p-{x) 



dx(p-{x) 



The fact that the curvature part is cubic in densities is 
related to the energy of the Fermi sea for a quadratic 
spectrum in ID. Moreover, such cubic terms lead to a vi- 
olation of the particle-hole symmetry, Lp± — > — </?±, again 



a consequence of the non-linearity of the spectrum away 
from half-filling. Finally, notice that the curvature term 
in Eg. jfTTI) agrees with the one found in the literature 
[lilllfTsee Appendix El for conventions on notations. 
In terms of the total phase and momentum this yields: 



F =1 j dx[nl + [d^^f] 

f dxd,^ [3Ul + 
bm J ^ 



{dxV? 



(18) 



According to phenomenological bosonization the added 
terms, {d^'p^)'^ and d^i^ 11^, are the lowest-order particle- 
hole symmetry violating terms compatible with both the 
U{1) invariance of the Hamiltonian, i.e. ~^ Lp -\- a 
where a generates the translations, and the space-time 
symmetries of the problem: is an odd function of x and 
an even function of t: 

P^--i=dx'P, j = dt(p, (19) 

where p is the total density (the smooth part of the den- 
sity close to (7 = 0; the 2kF part has been neglected), 
j the total current of the fermions and 11;^ ^ dt<p, see 
Eq. p2|) below. Notice also, from Eq. p9| . that the orig- 
inal (spinless) fermions correspond to kinks in ip. 

To include forward scattering interactions we follow 
the usual prescription of re-scaling the fields: 



7p V- 



ip, V 



in the Tomonaga-Luttinger part of the Hamiltonian, 
where: 



7p 



'l + y4,c/2-y2,c/2 



y l + y4,c/2 + y2,c/2' 
Up = v^{l+y4,c/2f-{y2,c/2)^, 

and Tji^c = 9i,c/''^v is the dimensionless coupling constant 
in the charge sector. In the following we take: 2/2. c = 
2/4, c — Vc which yields: 



1 



Ip 



V \/l + yc 



(20) 



The interactions affect the free fermion Hamiltonian den- 
sity as: 

1 

— (a^ipr - „ 

7p 6m 



n 



7p h2 + - {d^^y 



d.p> [ml + {d^ipf] 



(21) 

To proceed further, we find it more convenient to go to 
the Lagrangian representation, —^8x9, which yields 
the following Euclidean action: 



Se = drdx 



Uplp 



6 



dxP] {dxOf 



-idxddrp + [dx'pf 
^Ip 



m 



7 {d^vf 



5 



where m! = Qm/^/n and r = it is the imaginary time. 
Integrating over the dual phase field 6 yields: 



1 



1 



The non-trivial denominator comes from the conjugated 
momentum which reads: 



Upjp (1 



dx^) 



(22) 



Expanding the previous expression in the lowest mean- 
ingful order in 1/m yields: 

2v ^ 

+ [3 {drvY - {dxvf] , (23) 

mv 

where we have used the fact that: 7pUp = v, independent 
of interactions22,. 

As an alternative to Eq. (|23[) . we may work with the 
complete action: 



2 27p 



m m 



(24) 



which does not require any expansion in 1/m. 

With the help of Eq. or Eq. (gU we wish to de- 
termine the effect of the irrelevant terms appearing due 
to curvature on the correlation functions of the model. 
In particular the long-distance part (close to g = 0; 
the 2fci?— part is neglected) of the polarization operator 
reads: 

n(a;,r) = - ( dM^.r) 9,^(0,0) ), 

TT 

where we have taken into account the spin degeneracy. 
From its expression in Fourier space all we need to com- 
pute is the phase correlator V: 

2q2 

n{iuj, q) = V{iuj, q), V{iuj, q) = {\(p{iuj, q)\'^). 

TT 

(25) 

Curvature terms make the actions non-Gaussian a sig- 
nature of the fact that they induce interactions between 
the bosonic excitations. We will treat them in pertur- 
bation theory in 1/m. Notice that for Eq. ((23)) curva- 
ture generates an infinite number of irrelevant terms and 
we have taken into account only the lowest order terms. 
This means that our treatment of curvature is essentially 
perturbative (even if we manage to re-sum all terms gen- 
erated by Eq. ([231) )• At second-order, both models of 
Eq. (|23p and Eq. ([M)) are equivalent as will be shown in 
more details below. 



B. Dynamical structure factor 

We now compute Vliujjq) via bosonization, i.e. in 
perturbation theory in 1/m. From Eq. (|23p the zeroth 
order reads: 



V<^°\iu,q) 



Uplp 



{iujy ~ {upqY ' 



from which: 



TTUp {iuY - {upqY 



(26) 



(27) 



which agrees with the zero order term of Eq. (jlOp in the 
free-fermion limit: 7p = 1 and Up — v. 

In the Euclidean action the cubic perturbation reads: 



__i 1_ f dqidq2dq3 

m'v^ 13^ J (27r)3 



S[l + 2 + 3] 



[3 giW2W3 - 9i'72'73] <Pc^i (^Ov'^s (92)¥'^3 (^s)- 

As will be discussed in the following the knowledge of 
the self-energy associated with this perturbation (at least 
in the lowest order) does not imply any meaningful re- 
summation of band-curvature effects at the level of the 
Dyson equation. We will therefore simply calculate cor- 
rections to the bare Bose propagator order by order: 



V{iuj, q) = X>(°) {iu, q) -f P^^) (j^^ ^) ^ 



(28) 



and similarly for the polarization operator. The lowest 
order correction is of the second order in band-curvature. 
It reads: 

V^''\iLO,q)^V^''\iu,q) ^{iuo.q) V<^'^\tto , q) , 

where E(itij, q) is the second-order self-energy part and 
is used without any re-summation implied. From the 
expression of the free bosonic Green function, Eq. (j26p . 
it will allow us to compute the second-order correction 
to the polarization operator, Eq. (^5)) : 



2-7^ 

U^^){tco,q) = -=^ 



{Upqf 



TT [{iujY - (upq) 



212 



^iiLU,q). (29) 



For the iy9— Lagrangian the diagrammatic theory assigns 
a solid line for the bare propagator of Eq. and the 
second-order self-energy part is displayed on Fig. [3l To 
compute this self-energy part one has to be careful in 
Wick ordering the 8-point correlation functions generated 
by the cubic action at second order. This care is due to 
the fact that, in the 1^9— representation, the vertices carry 
products of frequency and momenta in a non-symmetric 
way which complicates the calculation of the combina- 
torial factor. The latter is crucial as we will see in the 
following. 

As an alternative to the 95— representation we may 
work in the ip, 0— representation. As we have noticed 
m Sec. Hira the ip-La grangian of Eq. (|23p is derived 
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FIG. 3; Second order diagram of the bosonic field theory for 
the Lagrangian of Eq. (|23p for spinless fermions. 



by neglecting an infinite number of irrelevant perturba- 
tions. On the other hand no such truncation appears 
in the 6*— Lagrangian of Eq. ([M)) . At the second or- 
der of the perturbation theory in curvature both ap- 
proaches should yield the same results while at higher 
orders the ip, 0— Lagrangian has to be used. Moreover, 
in the ip, 61— representation the vertices carry products of 
momentum in a symmetric way which simplifies the Wick 
ordering at the expense of having more diagrams to eval- 
uate. The latter can be seen from the Fourier transform 
of Eq. ((24| which yields the following perturbations: 



{^{iio', q)0{—iuj, —q)) 
{B{luj,q)ip{-iu,-q)) 



FIG. 4: The quantum field theory associated with the 
if, 6^— Lagrangian of Eq. (|24|l for spinless fermions. 




Se 



i 1 
m' 



E 



dqidq2dq3 
(27r)3 



(5[l + 2 + 3] 



The quadratic part of Eq. (1241) yields the zero-order cor- 
relators: 





= (|^(^c.,g)p)=- 








^{mc.,q)f)^~ 


Up 1 




= ('^(iw, g)6'(-iw, 


qy - 

q[{iuj'f - {upq)^]' 






(30) 



Notice that while the charge phase (p is related to 
particle- hole pairings, ~ Xl'"', its dual, 9, is a "su- 
perconducting" phase and v'j^g is related to particle- 
particle pairings. The fact that <p and 9 appear on equal 
footing in Eq. (p3)) is related to the absence of long-range 
order (either of the charge-density wave type or the su- 
perconducting one) in ID. 

The diagrammatic theory then assigns a solid line for 
the field and a double line for the 6*— field, see Fig. H) 

Expanding the perturbation up to second order and 
Wick ordering the 8— point correlation functions yields 
seven different types of diagrams represented on Fig. [S] 
They are classified according to the number of 0— wings 
attached to the self-energy parts. These diagrams, Fig.O 
yield exactly the same self-energy contribution as the one 
of the (p— representation. Fig. [3l 

The final result for the small q self-energy reads (at 



SLBl 



qiq2q3 [3Vu,Aqi)0u;:,{q2)9u;:,{q3) + V'^i(gi)¥'(^2(<72)¥'(^3('73)] 



SLCl 



FIG. 5: Second order diagrams of the bosonic field theory 
for the if, 61— Lagrangian of Eq. (|24|l for spinless fermions. SL 
stands for: Spin Less. 



T = for simplicity): 

^^pi^,q) = 

3 " 



96uo 



+ 



7> 
64m. 



1 + — 
1 

7p 



(31) 



— ) Upq [S[lu - Upq] - S[uj + Upq] 



[u? - {upqf] Tp{uJ,q), 



where we have restored m' = Gm/y^ and defined the 
kinematic factor: 



Tp{uJ,q) = 0[uj - Upq] - 9[-uj - Upq]. 



Multi-pair continuum 



(32) 



Away from the plasmon-cone, uj » Up\q\, Eq. ((51]) 
yields the small wave- vector, high-frequency part of the 
DSF: 



S^^H^,q) 



2\2 



(l-7p^ ) 

32Un 



2\ 2 

m 



— \ — Tp{uj,q), 

LO^ - {UpqY 



(33) 



7 




FIG. 6: Fourth order diagrams of the bosonic field theory for 
the (/p— Lagrangian of spinless fermions. 



with the help of the second order self-energy. Indeed, 
such a partial re-summation does not cure the i5— function 
singularities at the new frequencies of Eq. Higher 
order corrections, some of which are represented on Fig. 
[6l have to be taken into account in order to develop a 
more efficient re-summation scheme. This important is- 
sue is however beyond the scope of the present paper. 



which is the leading contribution to the DSF in this fre- 
quency regime. This tail is by now well known in the 
literature, cf. Refs. I^l2lll23 |. and describes the decay 
of a plasmon into two plasmon excitations as can be seen 
from the vertices of Fig. H) In fermionic language the 
plasmon corresponds to a coherent particle-hole pair. As 
a consequence, the tail is termed the two-pair excitation 
continuum, see Figs. [1] and [21 Notice that the kinematic 
factor appearing in Eq. see also Eq. ([5^ . is due to 
the fact that, in ID, low-energy particle-hole pairs cannot 
exist for arbitrarily small momentum. 



Z. Single-patr continuum 

Close to the plasmon-cone, oj = Up\q\, the total polar- 
ization operator may be again conveniently re-casted in 
the Matsubara form: 



n{iuj, q) = n(°) {iu, q) + n(2) (iuj, q) 



(34) 



where the dots refer to higher order band-curvature cor- 
rections, c/. the fourth order contributions of Fig. [6l and 
from Eqs. and 



n(2)(zc.,g) = 



{u,qfr 
(35) 

We now perform a re-summation of band-curvature cor- 
rections from the second-order self-energy part. For 
q <^ kp, this yields new poles for the dressed polarization 
operator: 



'^±iq) = Up\q\ ± 



2 3_ 

3 4 



7p 



Ip 2m 



(36) 



We therefore recover the 2-parametric structure of the 
single-pair continuum. Moreover, up to a numerical pref- 
actor of the order of unity, the width: 



5uj cx , 

TO* 



TO, 



(37) 



increases with increasing repulsive interactions. This 
agrees qualitatively with the exact results of Refs. [l^. 

On the other hand, the precise line-shape of the DSF in 
the single-pair continuum (divergency on the lower edge 
and crossover to the multi-pair continuum at the upper 
edge, see Fig. [5]) cannot be accessed, even qualitatively. 



IV. SPINFUL FERMIONS 

A. Model and spin-charge coupling 

Following the previous sections we now move to spinful 
fermions by introducing an additional internal, spin-1/2, 
degree of freedom to the fields: cr =t) ii so that: 



1 



'2ti 



exp 



2 l^^T 



X 

dx' Haix') 



This yields: 



_2y/n 
3to 

with the densities: 



dx [{d^^+^^f + (5,<P_,,)3] , 



We then introduce the charge and spin fields together 
with their canonically conjugated fields: 



1 



1 



Following Ref. J^], the Hamiltonian density then takes 
the following form: 



n 



+ 



7pn^ + - {d^ipY 

7<T 



nc 7 [{d^^f + 3 d^ifi u 



to 



(38a) 
(38b) 

(38c) 



+ 3 d.ifi [id^af + nl] + 6 9,(7 n,] , 

where to' = 6m/ ^/tt. In Eqs. p8p interactions in the 
charge and spin sectors have been included. For the latter 
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we have introduced the following parameters: 



la 



/ l + y4,s/2~2/2,s/2 

Y l + y4,s/2 + 2/2../2' 
«y^(l+y4,./2)2-(y2,./2)2 



where yi^s = gi,s/''^v is the dimensionless coupling con- 
stant in the spin sector. Just as for the charge sector in 
the following we take: ?/2,s = 2/4, s — Us, which yields: 



1 



la = 



(= 1), 



= vVT+Vs (= «)• (39) 



In the absence of curvature (m' oo) only Eq. (j38b|) 
remains which corresponds to the usual Tomonaga- 
Luttinger model. It's main feature is the separation of 
the Hilbert spaces associated with charge and spin sec- 
tors (the spin and charge parts commute with each-other; 
recall the anomalous commutations: [(^^(a::), n(p(x')] = 
iS{x — x'), and similarly for a and Ho-, while other com- 
mutators are zero). Moreover, as can be seen from: 

Ppix) = ^\f~ ^^^^ Pai^) = ^\[~ (40) 

- dtcr, 

TT 



jpix) ^ \l - dt(p, ja{x) 



kinks in if determine the density of fermions while, inde- 
pendently, kinks in a determine the magnetization den- 
sity due to the spin of these fermions (we are consider- 
ing the smooth part of the densities close to g = 0; the 
q — 2kF parts are neglected). Charge and spin travel 
at different velocities: Up and u^-, respectively, which are 
determined by the interactions among the fermions. As 
a result, charge flies away from spin so that the interact- 
ing fermions decompose into two elementary excitations. 
These are the spinless charged holon and anti-holon and 
the chargeless spin-1/2 spinon, see the textbooks 
Such elementary excitations do not appear explicitly. 
Only pairs of these do. They correspond to the physi- 
cal excitations of the system, see also [oHTol . In the spin 
sector they correspond to spin-singlet and -triplet excita- 
tions which are degenerate (they consist of two chargeless 
spin-1/2 spinous). In the charge sector they correspond 
to the charge-neutral holon- antiholon excitations other- 
wise known as particle-hole pairs. The Gaussian nature 
of the TL model implies that these excitations are free. 

In the presence of curvature the physical excitations ac- 
quire a finite lifetime. In particular, the first two terms in 
Eq. (|38cp are identical to the ones found in the spinless 
case. They lead to the decay of the charge excitations 
(particle-hole pairs) within the charge sector. More im- 
portantly, the next terms in Eq. (j38cp show that spin 
and charge Hilbert spaces are no more decoupled when 
curvature is introduced, see Ref. [ij]. This opens new 
channels for the decay of the physical excitations. In- 
deed, the third and fourth terms couple the charge den- 
sity to the spin density and current, respectively. And 



the last term describes the coupling between the charge 
and spin currents. Moreover, we see from Eq. pScp that 
spin and charge degrees of freedom do not enter in a 
symmetric way. As a matter of fact, the decay of spin 
excitations can only take place by affecting the whole 
charge density along the wire. These field theory argu- 
ments are in qualitative agreement with results obtained 
from more microscopic approaches. As we know from the 
exact solution of the Hubbard model, e.g. Refs. [olfToj for 
reviews, the spin-charge decoupling exists at the level of 
the elementary excitations^!: holons and spinous. How- 
ever, below half-filling, physical spin excitations involve 
a rearrangement of charge degrees of freedom and vice 
versa^^. Only in the low-energy limit, i.e. in the limit 
of a linear single-particle spectrum, do spin and charge 
sectors decouple completely. 



Even though we have derived Eq. p8|) the general form 
of the Hamiltonian and the particular asymmetry be- 
tween spin and charge may be understood from general 
symmetry considerations. This allows a phenomenologi- 
cal approach to the model. As in the spinless case, the 
fact that all curvature terms are cubic is related to the 
energy of the Fermi sea for a quadratic single-particle 
spectrum. It implies that the particle-hole symmetry is 
violated which is natural away from half-filling for such 
a spectrum. Moreover, the curvature terms in Eq. (j38cp . 
are compatible with the basic symmetries of the prob- 
lem: (p and cr are odd functions of x and even functions 
of t, see Eqs. pO|) . Finally, Eq. (I38cp is invariant un- 
der U{1) transformations of the charge fields and SU{2) 
transformations of the spin fields. The latter is not quite 
explicit here since we are using the Abelian bosoniza- 
tion technique^'^. Nevertheless, as a weak manifestation 
of the SU {2)i symmetry, all spin-fields in Eq. (|38cp are 
paired, e.g. a term of the form dx^p^ip^a is compatible 
with the basic symmetries but is not SU{2) invariant. 
Another consequence of the SU (2) spin-rotational invari- 
ance of the theory, e.g. see Q, is that ys = —yi/2 where 
yi is the dimensionless coupling constant arising from 
backscattcring. Thus, from Eq. (jSQ]) and in the absence 
of backscattering^^, the spin parameters are the free ones: 
Ua = V and 7o- = 1. We will work with general values of 
these parameters in what follows and set them to their 
non-interacting value at the end of the calculations. 



B. Field Theory 



To proceed further, we go to the Lagrangian formula- 
tion by introducing the dual fields: 9p and such that: 
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[^{iuj,q)\'^) 

6p{iuj, q)ip{—iw, -q)) 

a(iuj, q)B„{-iw, -q)) 
d„{iuj,q)a{-iuj, -q)) 




Tjqiq'iq-i 



7919293 



FIG. 7: The quantum field theory associated with the La- 
grangian of Eq. pif) for spinful fermions. 



Hip = dxOp and — dxOa- This yields: 
2 Zjp 



1 

TO 

3 

TO 

6 



- dx<p [{dxvf + i{dxepf 

- [{dxcrf + {dxO^f] 



dxOp dxCT dxOa 



(41) 



The field theory associated with Eq. pT|) is defined on 
Figlll In the charge sector, bare correlators are given by 
Eqs. ([30|) with 9 = 9p. Following the spinless case the 
correlator is associated with the particle-hole pair 
excitations. These are charge-neutral excitations, even 
though we will often refer to them as charge excitations, 
for simplicity, as holons and anti-holons never appear ex- 
plicitly. Moreover, as there is no spin-charge coupling at 
the level of the bare propagators, these charge-neutral 
excitations are spinless. In the spin sector, the bare cor- 
relators read: 



I?W(^c.,g) = (|a(^c.,g)p) 



v'■^U^u^,q)^{\9,{^uJ,q)\^)=-^ 



7cr - {uaqY ' 



= {(^{i^,q)9a{~i^,-q)) 
Vll{^^,q)^v(^l{^u:,q), 



{Hiuj,qr 



FIG. 8: The quantum field theory related to the Lagrangian 
of Eq. (|43p for spinful fermions. 



Alternatively, we may also integrate over the dual fields 
in Eq. ([iT|) and work with the following reduced La- 
grangian: 



J^E [<^, ct] 



1 

2> 
1 

2^ 
1 



1 



(drip)^ + Up [dx^f 



— {drcrf' + u„ idx<j) 



dxip [3(9.(^)2 - v^dx^)^] 
6 

H r-j dr^ dx(J dra, (43) 



where only the lowest order terms in 1/to' have been 
taken into account. The corresponding field theory is 
given on Fig. [8] 

With the help of the models of Eqs. (gl]) and (gSl), we 
will determine the effects of the spin-charge mixing terms 
arising from band-curvature on the response functions of 
the system at small wave- vector q. There are two of 
them: the charge DSF and the spin (or magnetic) DSF 
(mixed correlators of the type: {if{iu}, q)a{—iuj, — g)), are 
zero within the present field theory because there is no 
spin-charge mixing at the level of the bare propagators). 



C. Dynamical charge structure factor 

We first compute the second-order corrections in band- 
curvature to the charge DSF of spinful fermions. As for 
the case of spinless fermions, it corresponds to the imag- 
inary part of the related correction for the polarization 
operator: 



n(2)(zc.,g) = -^- 



{upqf 



(upq) 



212 



Ep(^c^,9), (44) 



g[(iw)2 - (u<^g)2]' 
(42) 



where P^ct is associated with spin-singlet or triplet exci- 
tations. We shall often refer to them as spin excitations, 
for simplicity, as spinous never appear explicitly. They 
are chargeless. 



where the self-energy part corresponds to the diagrams 
of Fig.[9]in the </3, cr— representation. The latter include a 
contribution which is similar to the spinless case, the first 
diagram in Fig. [9l arising from the decay of a charge (- 
neutral) boson into two charge (-neutral) bosons, see the 
first vertex in Fig. [51 Because of band curvature, spin- 
charge separation is violated and spin degrees of freedom 
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FIG. 9: Second-order diagrams contributing to the 
charge density correlation function generated from the 
ifi, cr— Lagrangian of Eq. (|43p for spinful fermions. 



1^ 



1^ 



and u = v), the diagrams of Fig. [9] reduce to the one of 
Fig. [3] with an additional factor of 2, the spin degener- 
acy. Similarly, the diagrams of Fig. [10] add with those 
of Fig. [5] and yield an overall factor of 2. To appreciate 
the spin-charge coupling one has therefore to impose an 
inequivalence between these degrees of freedom. 

This arises naturally for arbitrary values of the param- 
eters, 7p ^ 7ct. The total self-energy, close to g = 0, of 
interacting spinful fermions then reads (at T = 0): 



p 

„2 



75 



(45) 



7> 



3 ■ 


2 


^. 




1 






+ 


1 " 


7 


2 

P J 


1 " 





Upq [5[uj - Upq] - 5[uj + Upq]] 



n 2 



7<t7p 

2 



where the kinematic factor reads: 



Uaq - Uaq] - + u^q]] 
\mup J 



\-uj - u„q\. 



(46) 



FIG. 10: Second-order diagrams mixing spin and charge 
degrees of freedom and contributing to the charge density 
correlation function of spinful fermions. These diagrams are 
generated with the </p, cr, Sp, So-— Lagrangian of Eq. (|4ip . The 
self-energy parts only depend on the spin degrees of freedom. 
CD stands for Charge Density. The set CDA is identical to 
the diagrams appearing in the spinless case, Fig. [5] and is 
not displayed. Notice that, up to the fact that the self-energy 
depends only on the spin degrees of freedom, these diagrams 
are the same as those of the spinless case. Fig. (5] 



provide another channel for the decay of the charge bo- 
son. This channel corresponds to the second vertex of 
Fig. [8] whereby the charge boson decays into two spin 
bosons (each corresponding to a spin-singlet excitation) . 
The corresponding, second-order in band-curvature, dia- 
gram is the second diagram of Fig. [51 

In the (/?, cr, 6*^— representation of Eq. (|4T|l the new 
diagram is given by the diagrams of Fig. [TUj where the 
same classification as in the spinless case. Fig. [Sj has 
been used, i.e. diagrams are classified according to the 
number of 0— wings attached to their self-energy parts. 
The increased number of diagrams in Fig. [TOj with respect 
to Fig. [9] reflects the non-trivial nature of the vertices of 
the cubic field theory under consideration. 

The diagrams of Figs. [Sj and [TUj have the same struc- 
ture as those of the corresponding spinless case, Figs. [3| 
and [H respectively. In the case where charge and spin 
parameters are equal (7p = 7cr = 7 and Up — ~ u)^ 
the free-fermion case being a special sub-case (7 = 1 



This general result shows that the density correlation 
function of spinful fermions has two peaks: a charge peak 
at = Upq but also a spin peak at w = u^q so that 
part of the charge spectral weight is carried by the spin- 
singlet. The additional spin peak to the DSF is a witness 
of the spin-charge coupling due to band-curvature. It 
arises from the opening of a new channel for the decay of 
charge (-neutral) excitations, i.e. into two gapless spin- 
singlets traveling at < Up. The sharpness of this spin 
peak translates the coherence of the rearrangement of the 
background charge due to the spin-singlet. 

In the absence of backscattering, i. e. for 70- = 1 and 
Uct = V, Eq. (|45p simplifies as: 



7p7r 
96u„ 



TT 

24^ 

7p7r 
"64u„ 



3 

1 + - 

1 ' 

1 + — 

7p. 



7p^ 



— Upq [S[lu - 
m ' 



(47) 

5[uj + Upq]] 



— vq [5[uj — vq\ — 5[ijJ + wg]] 



[uj''~{upqf] Tp{LU,q). 



1. Multi-pair continuum 

Eq. 1^5)) shows that both charge and spin peaks have 
tails due to the interactions between the fermions com- 
posing the bosonic excitations. For ld ^ Up\q\ and from 
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vq 



Upq 



FIG. 11: Schematic view on the low-momentum [q <C kp) 
part of the spectrum of excitation of spinful fermions and 
the corresponding charge DSF. Interactions are present in the 
charge sector (7p < 1) while the spin sector is free (70- = 1)- 
Because of the violation of spin-charge separation due to 
band-curvature two peaks appear. The higher one corre- 
sponds to the usual charge peak and has a high-frequency 
tail (displayed in light gray). The lower one is the spin peak. 
Within bosonization, the width of the spin peak and the pre- 
cise line-shape of the DSF are not accessible. 



Eqs. (|45p and (|44p . the leading contribution to the charge 
DSF reads: 



5f (-,9) = 



(48) 



32up 

(l-7g)^ 
2,2u„ 



q 

m J Lu 
m 



Notice that, in the frequency range: Ua-\q\ <C ijJ ^ Up\q\, a 
small spectral weight is provided by the spin excitations. 

In the limit where backscattering is neglected, Eq. (|47)) . 
the spin peak at lo = v\q\ has no long-range tails and for 
strong repulsive interactions, 7p <C 1, is well separated 
from the charge peak at a; = Up\q\ = v\q\/jp ^ v\q\. 



In Eq. ([50)) the second term arises from the spin charge 
coupling. If we perform a re-summation of the second- 
order results, the polarization operator acquires new 
poles at the frequencies solving: 



(51) 



37p 



3 + 7' 



This equation is second order in u)"^ — {upqY from which 
we recover the 2— parametric nature of the charge peak. 
For q/kp ^ 1, the latter reads: 



^{q) ~ K\q\ ± 



2m* 



-Oiiq/kpf), 



(52) 



where m* is still given by Eq. ([57)1 . At this order the 
presence of the spin peak does not affect the width of the 
charge peak. The main effect of the spin-charge coupling, 
within the present re-summation based on the second 
order self-energy part, is then to shift the velocity of the 
charge excitations: u* ~ Up[l — (q/kp)^]. 



On the other hand, Eq. ([5T)) is first order in i 



{vqf 



This implies that, within a re-summation based on the 
second order self-energy part, the spin peak remains 
single-parametric with a velocity slightly renormalized 
by the band-curvature: v* « v[l — [q/kpY]- 

Our results are schematically displayed on Fig. [11] 
Within the present, perturbative, approach we cannot 
access the precise line-shape of the charge DSF in the 
vicinity of w w "pi?! ^-nd uj ~ v\q\. 



2. Single-pair continuum 

We now focus on the vicinity of the charge peak, 
LU « Up\q\ and neglect backscattering so that charge and 
spin peaks are well separated from eachother. Taking 
into account the leading contribution, Eq. ([77)1 . the to- 
tal polarization operator up to second order in band- 
curvature reads: 



n{iLU, q) = n(") {iuj, q) + n(2) {iuj, q) 



(49) 



where the dots refer to higher order corrections and, from 
Eqs. ([T7| and ([^1| . the second order correction to the 
polarization operator reads: 



DTT 



(50) 



2e. 

Au 



1 + - 

7p 



1 

7p, 



1 2 



{vqf{upqf 



[(iuY - {vqY][{iuY - (upq) 
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D. Dynamical spin structure factor 

We now focus on the spin degrees of freedom. The 
long-distance part (close to g = 0; the 2fci?— part is ne- 
glected) of the spin-density correlation function is defined 
as: 

X{x, r) = - {dMx, T)dMO, 0)) , (53) 

TT 

and, in Fourier space, reads: 



9) = 'D^„{iuj,q). 

TT 



(54) 



The dynamic spin structure factor corresponds to the 
dissipative part of this susceptibility: 

S,iLO,q) = -^x''{LU,q). (55) 

In the absence of curvature the result is well known: 

S"Aco,q)^^ {u^qf 6[co-u.\q\], (56) 

Urr 
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FIG. 12: Second-order diagrams contributing to the 
spin-density correlation function and generated from the 
</3, CT— Lagrangian of Eq. (|43p for spinful fermions. 



with obvious similarity with the dynamical charge struc- 
ture factor. 

Including curvature, the second-order correction to the 
spin density correlation function reads: 



where the self-energy is given by the diagrams of 
Fig. [12] in the the ip, cr— representation of Lagrangian 
Eq. ([43|l . These diagrams are built from the third ver- 
tex displayed on Fig. [8] which shows that a spin excitation 
(singlet, triplet) can decay only into a mixed spin (singlet, 
triplet, resp.) — charge-neutral excitation. Contrary 
to charge fluctuations which had two-channels through 
which they could decay (one of them not affecting the 
spin degrees of freedom and therefore surviving in the 
spinless case) the decay of spin fluctuations cannot pro- 
ceed without affecting the whole charge background. 



SDC3 




FIG. 13: Second-order diagrams contributing to the 
spin-density correlation function and generated from the 
</5, cr, Sp, Lagrangian of Eq. (|41|> for spinful fermions. All 
the self-energy parts mix spin and charge. SD stands for Spin 
Density. In the limit where spin and charge degrees of free- 
dom are equal these diagrams reduce to those of Fig. [S] up to 
spin degeneracy. 



In the (p, CT, ^p, representation of Lagrangian 
Eq. (|4ip , which reveals the non-trivial nature of the ver- 
tices of the cubic field theory, these mixed spin-charge 
self-energy parts are displayed on Fig. [T31 



In the case where spin and charge degrees of freedom 
are equivalent {jp — — 1 and Up = Ua- — u), the free- 
fermion case being a special sub-case (7 = 1 and u = v), 
spin and charge bare propagators are equal so that there 
is no more difference between the solid and dashed lines 
in Figs. [13 and [T31 One can then easily check that in 
this case the diagrams of Figs. and [THl are identical to 
those of Figs. [3] and [5l respectively, with an additional 
factor of 2 due to the spin degeneracy. 



All non-trivial effects that we shall observe below are 
therefore related to the fact that spin and charge excita- 
tions travel at different velocities {-jp ^ 70-). Even though 
calculations may be carried out for arbitrary spin and 
charge parameters we will focus here, for simplicity, on 
the case where backscattering is neglected so that the 
spin sector is free (7^ = 1 and = v). The self-energy 
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then reads (at T 



for simplicity): 



(58) 
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Notice that in the hmit where 7p ^ 1 in Eq. ([58|) the 
first two terms become (5— functions centered around ±vq 
and the last two terms vanish. On the other hand for 
7p < 1, the combined effect of repulsive interactions and 
second-order curvature effects broaden the peaks at ±vq 
by transferring spectral weight to frequencies reaching 
the charge frequency (up<?) with additional long-range 
tails. This is more conveniently seen on the magnetic 
DSF. The latter is non-zero all the way between the spin 
and charge frequencies: 

Si'Huj,q)^2!i [e[-,u + Upq]-e[-,u + vq]] 

Av \ mv / 



[l/7p + l]^ LO-Upq [l/7p-l]2 u + Upq 



|l/7p-iP LO-vq |l/7p + l|3 uj-vq 



, (59) 



where, for simplicity, we have set g > and focused on 
LO > vq. 

At LU = vq this expression diverges but only alge- 
braically and not as a (5— function as for the zero-order 
term or in the absence of interactions. 

At UJ = Upq the magnetic DSF is finite and equals: 

^^'^(""^''^^-FhT^ (-)^ 

2v |l/7p-fl|'' \mv/ 

so that part of the magnetic spectral weight is carried by 
the charge excitations. 

At large frequencies, uj ^ Upq, tails appear and read: 



l/7p + 1 



q- \ 6[uj- Upq] 
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FIG. 14: The magnetic structure factor, Sa(uj,q), at a fixed 
q, as a function of the dimensionless frequency Co — u;/vq 
and for different interaction strengths in the charge sector 
{ja ~ 1) • This figure shows a transfer of the spectral weight to 
higher frequencies as the interaction strength increases. This 
transfer originates from the spin-charge coupling due to band- 
curvature. The spectral weight then extends from the spin 
excitation frequency, (Do- = 1, to the interaction-dependent 
charge excitation frequency, cOp = l/7p, with additional high 
frequency tails. The peak at (D = 1 corresponds to the free 
fermion result: 7p = 1 {loct = uop = 1). The second curve 
corresponds to the second-order curvature correction at jp = 
0.9 (Cop = 1.1). The third curve corresponds to the second- 
order curvature correction at 7p = 0.1 (Cop — 10). 



Notice that the second term in Eq. (|6ip yields a non- 
zero dissipation at infinite frequencies. Nevertheless, the 
spin-structure factor satisfies the sum-rule: 



duj UJ Sa-{u}, q) = 



V q 



(62) 



but with a re-normalized spin- velocity, v* , which reads 

7p' [l-7p]' 



1 



4 |l + 7pp 



V < V < Uo 



(63) 



In Eq. (inH) a large frequency cut-off a,t ep — mv'^ has 
been included and the final result was obtained after 
sending this cut-off to infinity. The re-normalization of 
the spin-velocity by interactions in the charge sector is 
precisely due to the last term in Eq. (pT|) . It is a signature 
that the physical excitations carrying the magnetic spec- 
tral weight, in the tails of the spin-structure factor, are 
spin (-singlet or -triplet) excitations dressed by charge- 
neutral ones. 

We therefore see that band-curvature softens the singu- 
larities of the spin structure factor already at the second 
order of perturbation theory. Such an effect arises from 
spin-charge mixing and is more dramatic than for the 
charge DSF. The latter was also affected by the viola- 
tion of the spin-charge separation due to band curvature 
(c/. the double peak structure, see Eqs. ([47]) and ((50|) and 
Fig. llip but its low frequency part was still (5— singular at 
second order in band-curvature. The sharpness of both 
charge and spin peaks in the charge DSF was due to the 
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decay of charge fluctuations into either a pair of charge- 
neutral excitations or a pair of spin-singlet excitations, 
both coherently traveling at velocities Up and Uc, respec- 
tively. On the other hand spin fluctuations decay via 
spin(-singlet or -triplet) — charge(-neutral) excitations. 
Both components of this pair lead to rearrangement of 
the background charge along the wire^i. But they travel 
at different velocities. As a consequence, the redistribu- 
tion of magnetic spectral weight is incoherent and affects 
an infinite number of modes. The magnetic DSF is plot- 
ted on Fig.[T4]for different values of the interactions in the 
charge sector. Because of the spin-charge coupling just 
described more spectral weight is transfered from the spin 
peak at Lo — vq to the charge peak at lu = UpQ = vq/^p, 
upon increasing interactions. 



V. CONCLUSION 

Using the bosonization technique we have studied the 
charge and spin equilibrium dynamics of ID fermions 
in clean quantum wires with forward scattering inter- 
actions, i.e. below half- filling, and band-curvature, i.e. 
quadratic, corrections to the linear single-particle spec- 
trum, cf. Eq. ([3]). The dynamics were accessed by cal- 
culating the charge and magnetic dynamical structure 
factors at small wave-vector q. Charge dynamics may be 
accessed experimentally with the help of energy and mo- 
mentum resolved spectroscopies, [231, electronic Raman 
spectroscopy, [2§|, or electron energy loss spectroscopies 

The introduction of band-curvature corrections al- 
lowed us to explore the high-energy physics of fermions 
in quantum wires. This a basic step beyond the stan- 
dard Tomonaga-Luttinger model which assumes a linear 
dispersion, Eq. ([T]). In solid-state physics this relativistic 
approximation is valid only at low-energies. It is also a 
basic step towards the Hubbard model (below half-filling) 
which has the full non-linear spectrum of lattice fermions, 
Eq. 0. 

Based on recent results for fermions without spin 
[20II2H we have discussed the efficiency and limitations 
of bosonization in presence of curvature corrections, cf. 
Sec. mil Our main result is Eq. ([3T|) for the self-energy 
part of the DSF, in second-order in band-curvature and 
to all orders in interactions. From this self-energy the 
2— pair continuum of the DSF has been derived, follow- 
ing Ref. [2ll|. and in agreement with other known results, 
cf. [16,22]. Moreover, a re-summation of band-curvature 
effects based on the second-order self-energy enabled us 
to emulate the 2-parametric structure of the single-pair 
continuum, cf. Eg. ()37p . in qualitative agreement with 
exact results [20||21 |. On the other hand the precise line 
shape of the DSF could not be derived with the help of 
this self-energy and a more efficient re-summation scheme 
has to be developed. 

In the case of fermions with spin we have focused on 
the spin-charge coupling which arises from the interplay 



between curvature and interactions, cf. Sec. IIVI follow- 
ing Ref. [l3| ■ The spin-charge coupling manifests by the 
opening of new channels for the decay of the excitations. 
For charge excitations, besides the usual decay into two 
other charge-neutral excitations, as in the spinless case, 
a decay into two spin-singlets takes place. On the other 
hand there is a unique channel of decay of spin (-singlet 
or triplet) excitations: via a pair composed of the spin 
excitation itself and a charge-neutral excitation. The de- 
cay of spin excitations therefore always affects the charge 
degrees of freedom. 

For the charge dynamics our main result is Eq. (j45p 
the self-energy part of the charge DSF, in second-order 
in band-curvature and to all orders in interactions. This 
equation and the corresponding Fig. [Tl] show that the 
charge DSF has a double peak structure: a charge peak 
at LO — Up\q\ but also a spin peak at uj = v\q\. The 
charge peak corresponds to an incoherent continuum of 
single spinless charge-neutral excitations. It has a high- 
frequency tail corresponding to the incoherent continuum 
of pairs of these spinless charge-neutral excitations. On 
the other hand the spin peak corresponds to the inco- 
herent continuum of chargeless single spin-singlet excita- 
tions. A re-summation based on the second order self- 
energy has revealed that the 2-parametric structure of 
the charge peak is qualitatively similar to the case of 
spinless fermions. The width of the spin peak as well 
as the precise line-shape of the charge DSF could not be 
derived with the help of this self-energy. For the spin dy- 
namics our main result is Eq. (j58p . the self-energy part 
of the spin DSF, in second-order in band-curvature and 
to all orders in interactions. This self-energy displays a 
transfer of magnetic spectral weight to higher frequen- 
cies. This continuous transfer extends from the spin fre- 
quency LO = v\q\ to the charge frequency to — Up\q\ with 
additional tails at lo ^ Up\q\. Because, Up — v/jp is 
interaction-dependent the transfer of magnetic spectral 
weight to higher frequencies increases with increasing re- 
pulsive interactions, see Fig. [HI 
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APPENDIX A: NOTATIONS 

Our notations are different from the ones often found 
in the literature, cf [UlllH. They have: 



1 



127rm 



with: ipR(L)ix) f.±ii':FX-zipR^L)_ We have therefore the 
following correspondence between these notations and 
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our notations: 



(Al) 



We may check that, indeed, this yields the same numer- 
ical factor: (2-y7r/3m)(l/87rY^) — l/12nm. 

Moreover, the chiral basis, Lp± , is often used. We which 
to relate the general irrelevant operator expressed in the 
ip basis: 

with its expression in the chiral basis {n — 8/7^ in the 
text). In real time this can be done with the help of: 
{l/v)dtip = = -dx(p+ + dx(p- and dx(p = + 
dx^P-, and yields: 

dx(p [nidrff - v'^{dx'pf] = 

-(n + 1) [{dxf+f + [dx^-f] + 

+ {n - 3) [{dxP+)'^ dxW- + dxP+ [dx^-f] ■ 

The case n — 3 corresponds to free fermions: 



dxf [3(9r<p)^ - v'^{dx(py 



-4 [{dxP+f + {dxP-f] 

- [{dx'PRf + (dx'PL)'] , 



and is the sum of two chiral terms. This implies that tails 
are in the crossed chiral terms. They can be isolated by 
a special (unrealistic) choice of the interactions: n = —1, 
which yields: 

dxp[{drpf+V^{dxpf] = 

4 [{dx(p+)'^ dx(p^ +dx(p+ {dx(p-f] = 

— 7= [{dx'PR)'^ dxipL - dx(fR {dxPL)"^] , 
7rv47r 

where Eq. (jAl[) has been used. 

APPENDIX B: OPERATOR PRODUCT 
EXPANSIONS 

Bosonization rests on the following vertex operator: 
1 



2-K 



exp 



±iv47r ip±{x) 



relating the chiral fermionic fields to the chiral bosonic 
fields. One then needs to express various operators (den- 
sity, Hamiltonian...) made out of fermions in terms of 
the bosonic fields. The operators contain the products of 
fermionic operators at coinciding points and a regulariza- 
tion procedure is required in order to remove the short- 
distance divergences. This can be achieved by point- 
splitting the operators, e.g.: 

^l(x)ipr,{x) = lim^^o [i^lix + e)iJr,(x) - {i/jlix + e)ip^{x))] 



where rj = ± and: 



(4(a;-f e)V±(x)) =±^. 



(Bl) 



The product of vertex operators is normal ordered ac- 
cording to: 



where: 



{lP^{x)lp^{x)) ^ \n[i{x-x)]. 



(B2) 



Notice that the average values in Eqs. (jBip and (|B2[) were 
taken over the TL action, i.e. no curvature correction 
was taken into account. We believe this approximation 
is sufficient, within the present perturbative approach, to 
bosonize the theory. Indeed, the first curvature correc- 
tions to, e.g. Eq. (jB2[) . are in 1/w? while the fermionic 
Hamiltonian has a contribution in l/m. At the level of 
the current algebra this approximation implies that we 
keep the usual (ultra-local) commutation relations: 

[dxip{x),IV{x')] = idx5(x - x') + 

where the dots correspond to the neglected curvature cor- 
rections. 

Altogether this yields: 



Pnix) = tpl{x)iJr,{x) 
= sgn(77) hnie- 



The first OPE reads: 



27rie 



g-isgn(j))\/47re9i:ip^(a:) _ ^ 



(B3) 



More generally: 
sgn(?7) 



1 



2Tri{y — a;) 



from which: 



tPl{y)dxiPvi^) 



sgn(?7) 



1 



27ri(y — x)^ 



-isgn{ri)V4^[^p,,{y)-<p,,{x)] _ ^ 



1 



dxPnix) e~*^Sii(»?)V47r[^„(y)-y„(a;)]^ 



and: 



ipl(y)dxxipvix) 



sgn(?7) 



1 



TTi(y — x)^ 
2 



V^{y - x)^ 
1 

V^{y - 

-sgn(?7) 



y-x 



g-isg.n(ri)ViTT[Vr,(y)~Vr,(x)] _ 2 
dxPn{x) e^*''Sn('?)V47r[ip„(y)-Vr,(2:)] 
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The short-distance products arc obtained by setting: y = Expanding in e the diverging terms in 1/e and 1/e^ cancel 
X + e and expanding to the lowest meaningful order in e. out and the third OPE reads: 
The diverging terms in 1/e cancel out and as a result the 
second OPE reads: 



1 



It may be introduced into the Dirac Hamiltonian density: 



= v[: {d,^+ix)f : + : {d,^-{x)f :] , 



where we have neglected the terms which vanish upon 
integrating over space. 

Similarly, the curvature terms read: 



ipl{x + e)da:xipr,ix) = sgn{r]) —r^ 

g-isgn(77)\/47r[e9a!Vn(^)+«^/2 9ica,*'r,(a;)+e^/6 d^xxV-ni"')] _ \ 

2 



V'U^)^^V'r)(a;) = -isgn(?7) : (ax<y!?r,(a;))^ : a^^<^^(a;). ^t(a;)a^^^^(a;) = 1|5 : {d^ipr,{x)f : (B5) 

(B4) 1 

-3zsgn(?7) : dx^Pn{x)dxx^ri{x) : -i^l^ dxxxfr,{x). 



30F 



It may be introduced into the curvature part of the 
Hamiltonian density: 



1 
2m 
2v^ 



i'+{x)dxxtp+{x) + tpl{x)dxtp-{x) 
3^ [:{dx<p+{x)f: + :{dx<p-ix)f-] 



+ 



+ 



0re 



; dx<frj{x) e~*''Sn(';)V47r[eaxV,(a:)+£^/2 dxxVvi^)] 



+Sgn(r/) - {dx'Pr,{x)f e-*^Sn(r,)V4^ea.c,,(a;) 



where we have neglected the terms which vanish upon in- 
tegrating over space or by symmetry considerations [up is 

an odd function of x so that (dx'p)^ is even but dx^f dxxV 
is odd and does not contribute to the energy]. 



J 



APPENDIX C: CORRELATION FUNCTIONS IN THE SPINLESS FERMION CASE 
1. Second-order perturbation in curvature: (^—representation 
The second order correlation function reads: 
^l'^^*^'^^!^ =2K^^ ^, , , J ^^<5[l + 2 + 3]5[l +2 +3] 

CJi ,a;2 :<-i^3 i<-^2'^3 

where Hp = 3/7^ and may be conveniently split into three contributions: 



dq.dq^dq^ dq[dq',dq'^ S\l + 2 + 315M' + 2' + 3'1 

Wi,W2,W3,wi,W2,W3 



«-"l^>^'^5^S?tlPi^ 



(27r)3 (27r) 

9l92g39i«293 ('^u;(5)<^-w(-g)<Pwi(9l)V'a;2(92)v'a;3(93)</'a>;(Q''l)</?wi(92)</?t^^(93))0, 



Mi^.q)?)f 



1 



E 



dqidq2dq3 dq[dq!2dq'3 



qiu)2U)sq'i^Wd. (</fw(g)</f-w (-g)</fwi {qi)vuj2 {q2)'Pw3 (93)9?^ {q'i)fui'^ {(^2)^^'^ (93))o 



5[l + 2 + 3](5[1' + 2' + 3'] 



npWp 1 



(m'M2)2 ^6 



E 



J dq^ djA^ ^ 2 + 3]^[1' + 2' + 3'] 



gia;2W3«ig293 ('Pa;(«)</'-a;(-g)<i2a;i(«l)^W2(a2)¥'w3(53)¥'a;;(al)</^a^^(92)¥'a;^(93)>0- 
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The first contribution is easy to Wick order and yields: 

18^2^4 



[m Up) p ^ J 



The imaginary part of the corresponding retarded self-energy reads: 

Q / x2/2\2 



^Ai^^q) = — ( — ) i:^] Upq [S[uj - Upq] - S[uj + Upq]] 



16up \Up 



where J-"p is given by Eq. (|32p . The second correlation function is simplified as: 

^n^quj „ 1 .r-^ /" rfp 2 

[m'UpY p ^-^ J zn 



(rn'u^j^ p — ' J Ztt 

— ^ [X>(zw,g)] — p{p-q)v{v-uj)V{iv,p)V{iv-iuj,p-q). 

[m Up) p J ZTT 



All terms contribute to the on-shell part of the correlation function but only the first contributes to the off-shell part. 
The result for the self-energy part reads: 

3^2 /^\2 /2\2 

q) = — I — I Upq [5[uj - Upq] - 6[uj + Upq]] 

OUp \Up J \TTl J 

r,2 / \ 2 / \ 2 

: P I \ I 1 \ r 2 / n2 



16up \Up 



Lo'-iupqY Tp{oj,q). (C2) 



Finally, the third correlation function is simplified as: 

{Wu,q)\^)^l) =- ^,^";tf [^(*^'g)]' ^ E / ^p{p-qy{^~OJ)V{^u,p)V{^u-^u^^p-q) 

2AnpU^qui „ 1 f dp , 
( I 2\2 [^(^^'g)] flE / (P-g)(^-^) ^(«^'P)^(*^-*^'P-g)■ 



Again, only the first term (which carries the q^ factor) contributes to the off-shell part whereas both contribute to 
the on-shell part. The corresponding self-energy part reads: 

Q / \2/2\2 

ov«/. . „N iup ( v\ ( q ^ 



Q'Sc.(lj,(j) =— ( — ) ( — ) Upq [5[u:-Upq]-5[uj + Upq]] 

u^~{upqf] TpiLU,q). (C3) 



Aup \Up 

o / \ 2 / N 2 

dUp I V \ I q 



8up \Up 



Adding the above self-energy parts the factors (1 + Up)'^ = (1 + 3/jp)'^ and (3 — np)^ = 9(1 — l/7p)2 appear and the 
self-energy is given by Eq. (PT|) . 
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2. Second-order perturbation in curvature: Lp, 6'— representation 

Expanding the perturbation up to second order yields three contributions: 

UJ\ .UJ2 -I^S ^<^l -1^9 :'^3 

9192939^9293 ('/'c^(9)<^-c^(-9)<y5t^i(9l)<^t^2(92)</?c^3(93)'^c^;(9'l)'/'^(92)V^^(93)>0, 

,91 ^ f dqidq2dq3 dq[dq'2dq'3 ;cri i o , Qi;rrv , n' , on 



2m'2 /36 ^ y (27r)3 (27r)3 

cJi ,a;2 ,t^3 ,t^^ ,t^2 



9i92939i9293 (V'^(9)</'-t^(-9)</'a;i(9i)^'a;2(92)6'a;3 (93)9'^; (9i)^'(^^ (92)^*1^,^ (93))o, 

,31 f dqidq2dq3 dq^dq'^dq'^ ^^1,0, 01 , n> , on 
2. / ^rT^'^[l + 2 + 3]<5[l +2 +3] 



m'2 /36 ^ 7 (27r)3 (27r)3 

ci;i .UJ2 ,^^3 - t^i ,^^2 ''^a 

9i92939i9293 (<y5t^(9)<y5-t^(-9)</'a;i(9i)^'a;2(92)6'a;3(93)</'^;(9i)v'c^i(92)'^c^,^(93)>0- 

The first contribution is exactly the same as the A— term of the previous paragraph. 

Wick ordering the next contributions we classify the diagrams according to the number of 6*— wings, see Fig. [S] 
Diagrams with no 6'— wing are labeled as SLA (SpinLess A). SLAl corresponds to the previous A term. The next 
terms read: 

V^Il,{^c.,q) = ~18[V(^l{^co,q)r {^f i ^ | g /(^ _ ,)2 ^(0) (,,^ ^(0) _ ^ _ ^) 



Their self-energy parts reads: 



3^2 / g2 \ 2 

[uj^ - {upq)'^] Tp{uj,q), 



16up \m'up 

r,2 



'p 

9 / q 



(m' ) ^^^^ ~ "'''^^ ^ ^ "'"^^^ 

[w^ - [upqf] Tp{uj,q), 



2 



16up72 \m'up 



9 

8up \m'u 
Diagrams with one 6*— wing (SLB) read: 

V^^l,{^^,q)^72V^^^,{^u;,q)V^^^{^u;,q) (-^)' i J] p2(^ „ ^)2 ^(0) (^^^ ^(0) _ ^ _ ^) 

2'il,(-,9) = -72Pj;'(-,9)2?^''^(*-,9) (;^)' ^ E j ^ p'{p - q? T^^^li^-.p) - ^^,P - q)- 

The corresponding self-energy parts read: 



2up7, 



2 



( ^ I Wp9 [(^i^ - Upq] - 5[u! + Upg]] 
y I Up? - Upq] - 5[uj + Upq]] . 
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Diagrams with two 0— wings (SLC) read: 



V. 



(2) 



{iuj,q) 



v-iuj,p-q) 



The corresponding self-energy parts read: 



9 



4up7p 
9 

8wp7p 



— - Upq \5\uj — Upq] — 5\uj + u 
m! ) 

\ ^ 

TV Up) ^ r i r 



Upq [5[uj — Upq\ — 6[uj + Upq] 



The total contribution of each set read 



16u,, 



1 



1 

7p 

_3 

2w 



7p^ 



oq [S[uj - Upq] - 6[uj + Upq]] 



- (upqf] Tp{uj,q) 



2up7 



I'Up) L 

/ 2 \ 2 

^ \ ~] '"p'? [^k - Wp?] - (^it^ + Up?]] • 

p V"^ / 



7p 



(C4) 



It is interesting to see that all diagrams contribute to the on-shell self-energy; however, the final expression for the 
tail comes only from the SLA set, i.e. those diagrams with no 6'— wings. The set SLB with a single 6'— wing does 
not have any off-shell contribution whereas the set SLC with two 6*— wings has off-shell contributions which cancel 
each-other. Notice that if the latter had a net contribution to the tails, this contribution would have higher powers of 
frequency which would violate the f-sum rule. In dimensionless units these additional contributions are: cx lojupq for 
SLB and cx (lojupq)'^ for SLC due to the presence of a single 0— wing and two 0— wings. On-shell, i.e. for w ~ ±Upq, 
such contributions are ±1 which yield the correct signs in front of the delta functions. 
Adding all the above self-energies yields: 



1 



97; 



16up 



1 ■ 



7p 

q 



Upq [5[uj - Upq] - S[uj + Upq]] 



[w^ - [upqf] Tp{uj,q). 



(C5) 



Substituting m! — 6m/ yields Eq. pT|) which is precisely what has been found in the previous paragraph. 



APPENDIX D: CORRELATION FUNCTIONS IN THE SPINFUL FERMION CASE 

1. Charge structure factor 



With the help of the Lagrangian Eq. (|4T|) the second-order contributions to the charge phase correlator have three 
contributions: 



(|iy9(ici;, g)!^)!^' = Similar to the spinless 
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which is identical to what is found in the spinless case, 

















18 [V'^^l{lUJ,q)Y 


( ' 
Vm'/ 


) ^ 


E 


/■ rfp 

/ 27r 


2 / 

p vp- 


qf {iv, p) {lu-lUJ,p- q) 


36[p(")(za;,g)]2 


f- 

Vm'/ 




E 

1/ 


/ 2tt 


p^ip- 


if T^al i^^^ P) ^iel {i^-i^.P- q) 




f- 

Vm'/ 


f - 


E 

1/ 


[ dp 
1 2tt 


p^ip- 


q? 'Df^g^ {iv, P) V'^^Jg^ {iv -iu,p-q), 



M^LO,q)f)P = 

-nV^li^.,q)vf,^i^.,q) (5)'^E 
-nV^li^.,q)V%[^u,q) (5)'^E 



^ p^b - q? v^^]^{w,p) vf^^g^ {iv -iuj,p- q) 

dp 



2tt 



q \2 1 



m' /' /3 



27r 



-36 {^y -p^ jt p'^'-^^' 



iiy — — 17). 



Based on the similarity of the calculation with the spinless case we directly give the expression of the various self-energy 
parts: 



1 



16m, 



37| 
8m„ 



1 \ 2 / \ 2 



To 
3 



2 

/ 2 \ 2 

( — ) Uaq [S[uj - Uaq] - S[uj + Uaq]] 
\m' J 

Ucrq [S[uj - u^q] - S[uj + Ucrq]] , 



\m' 






1 " 


1 + 










(si 




V to' 





(Dl) 



where is given by Eq. (j46|) . Notice that these expressions satisfy the basic properties of the diagrams with zero, 
one or two 0— wings. 

Adding all the above self-energies and substituting m' — 6m/y/n yields: 



1 2 
1 + - + — 



7a TT 



7a 7a7p 

2 



Ua<7 [^[t^ - Ucrq] - 5[UJ + U^rq]] 



[uj^ ~ (u^q)^] T„{uj,q). 



(D2) 



Adding the part arising from the charge degrees of freedom, the total self-energy of interacting spinful fermions is 
given by Eq. (liS)) . 
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2. Spin structure factor 

With the help of the Lagrangian Eq. (|4T|) the second-order contributions to the self-energy part of the spin structure 
factor read: 



f- 




E 

u 


/•dp 

/ 27r 




g)2 V^°l{^,,,p) Vf}{iv-iu,p-q) 


f- 




E 


[dp 
1 27r 




if ^e°e^(«f''P) -iu},p- q) 


f- 

Vm'/ 




E 


/■dp 
/ 27r 







72 / <? \2 1 

7cr Ma? Vm'/ /5 

72 ILJ / q \ 2 1 

7cr Uo-q Vm'/ /5 

72 iuj / q \ 2 1 



7cr Uaq \m' J (3 
72 / q V 1 



7cr Wo-q Vm'/ /3 



^P^ip-qf vf^Q^ {iv, p) Vl^l {iiy~iuj,p- q) 
dp 



27r 

dp 
2^ 



^ P^{p-qf T^^^eji'^,?) 'Df^g^iiv -i(^,p- q), 



36 / iLo 



la \Uaq 

36 / icj 



(;^)' ^ E J ^p\p-qrm^-^p)^''X(^-~^-^p-i) 



72 f iuj \ / g \2 1 /• , • ^ 



il \u„q) \m' I P 



The corresponding imaginary parts read: 



lpl<y 



u„q 



\up - UaY'q^ 



u ~ u^q] [uJ - Upq] 
+ + Uaq] [t^ + Upq] 





Up[LU - Uaq]' 




'u„[uj- Upq]' 






- n.B 






Up-Ua 




Up-U„ 






Up[ui -1- Uaq] 




u„ [uj + Upq] 






- ri B 






Up-Ucr 




Up - Ua 





Ipla 



1 



1 



\m') 



m' J \up + Ua\^q^ 



[uj + Uaq] [(^ — Upq] 
+ [uo — Uaq] [oJ + Upq] 



riB 



riB 



Ua\OJ - U, 



.9] 



Up + U„ 

Ug- [lo + Upq] 
Up + u„ 



UB 



riB 



Up[uj + Uaq] 



Up + Ug 
Up[uj - Uaq] 
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where ub is the Bose occupation function. Summing these contributions yields: 
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At T = and in the hmit where spin and charge parameters are equal (7^ = 7p, m^- = Up) we recover the result 
of Eq. ([51]) with an additional factor of 2, i.e. the spin degeneracy. At T = and in the absence of backscattering 
(70- = 1, Ua = v), Eq. (jD6p reduces to Eq. 
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